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Abstract 

The deuteron electromagnetic form factors, A{Q'^) and B{Q'^), and the tensor 
polarization T2o(Q^), are unambiguously calculated within the front-form relativis- 
tic Hamiltonian dynamics, by using a novel current, built up from one-body terms, 
which fulfills Poincare, parity and time reversal covariance, together with Hermiticity 
and the continuity equation. A simultaneous description of the experimental data for 
the three deuteron form factors is achieved up to < OA{GeV/c)'^ . At higher mo- 
mentum transfer, different nucleon-nucleon interactions strongly affect A^Q"^), B{Q'^), 
and T2o((5^) and the effects of the interactions can be related to S-state kinetic energy 
in the deuteron. Different nucleon form factor models have huge effects on A{Q'^), 
smaller effects on B{Q^) and essentially none on T2o((5^)- 
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1 Introduction 



The deuteron is a fundamental system for our understanding of nuclear physics and a chal- 
lenge to our ability to describe nuclei as systems of interacting nucleons with a well defined 
internal structure, without an explicit use of their quark substructure. In particular elastic 
electron-deuteron scattering is a crucial test for deuteron models. 

There exists a wide literature (see, e.g., [|l], H, ^ H, ||, ||, |^ and Refs. quoted therein) 
devoted to the investigation of deuteron electromagnetic (em) properties and in particular 
to the accuracy of the one-body impulse approximation (lA) for the current operator. It 
is usually believed that effects beyond lA, as meson-exchange currents, A^iV-pair creation 
terms {Z graphs), isobar configurations in the deuteron wavefunction, etc. are important for 
the explanation of existing data. However, the contributions of these effects are essentially 
model dependent and, furthermore, obviously depend on the reference frame (see, e.g., 
Refs. [|, 0). 

Since precise measurements of the deuteron elastic form factors have been recently 
performed in a wide range of momentum transfer, up to = —q^ = 6(GeV/c)^ for ^(Q^) 
[111], O], theoretical models require a relativistic framework for a reliable description of the 



available data. Furthermore, it has been recently shown [|T3| that relativistic effects are 
relevant even for static deuteron properties, as the magnetic and quadrupole moments. 

An essential requirement for relativistic approaches is the covariance of the current 
operator with respect to Poincare group transformations. This requirement is non-trivial for 
systems of interacting particles, since some of the generators are interaction dependent. 

A widely adopted relativistic framework for the study of deuteron em properties is 
the front-form Hamiltonian dynamics (FFHD) with a finite number of particles (see Refs. 
T^ , |T5|] and Refs. [|l^, [T^ for extensive reviews), which gives the possibility to retain the 



large amount of successfull phenomenology developed within the nonrelativistic approaches. 
Indeed, in the FFHD seven, out of ten, Poincare generators are interaction free, in particular 
the boost generators, while P~ = (Pq ~ Pz)/V^ {P is the total momentum of the system) 
and the rotations around the x and y axes contain the dynamics. Only the two-nucleon state 
is usually considered and the wavefunction of the system factorizes for any front-form boost 
in an eigenfunction of the total momentum times an intrinsic wavefunction, depending only 
on internal variables. Therefore, in the case of elastic e — d scattering, one can express the 
three deuteron em form factors, determined by three independent matrix elements of the 
current, in terms of the deuteron internal wave function and the elastic em nucleon form 
factors (f.f.). 

In the FFHD the em properties of the deuteron were usually studied in the reference 
frame where g"*" = (qo+Qz)/ = (g is the momentum transfer) [jl], ^, |18], [1^, ^ The 
one-body approximation was used to define three matrix elements of the plus component of 
the current, while the other matrix elements of the plus component and the other components 
of the current were properly defined in order to fulfill Poincare covariance, Hermiticity and 
current conservation. However, for spin-one systems, as the deuteron, this procedure is not 
unique and gives rise to ambiguities in the calculation of the form factors 
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In Ref. [0, using a representation of the Poincare group within FFHD, we have 
shown that extended Poincare covariance (i.e., Poincare plus parity, P, and time reversal, T, 
covariance) is fulfilled by the current which has a one-body form in the Breit reference frame 
where the initial and final momenta of the system are directed along the spin quantization 
axis (gl = q — Qz^z = 0). Furthermore, we have shown that Hermiticity and current 
conservation can be easily implemented. An important feature of our approach is that it 
allows one to use the same definition for all the matrix elements of the current, without 
ambiguities. 

In a previous paper |T^, as a test of our current, we evaluated the deuteron form 
factors at = 0, namely the magnetic moment, fid, and the quadrupole moment, Qd, of 
the deuteron, which are not affected by the uncertainties in the knowledge of the neutron em 
form factors at finite momentum transfers. The deuteron magnetic and quadrupole moments 
represented a longstanding problem in nuclear physics. Indeed, theoretical calculations were 
not able to accurately reproduce in a coherent approach the experimental values for both 
quantities at the same time, although a variety of approaches have been attempted, by 
changing the tensor content of the nucleon-nucleon (A^ — A^) interaction, or considering two- 



body current contributions, both in non-relativistic and in relativistic frameworks [^4|, ^ ^ 



0. On the contrary, using our Poincare covariant current operator, this usual disagreement 
between theoretical and experimental results was reduced to 0.5% for fid and to 2% for Qd 
by using interactions able to reproduce the experimental value of the deuteron asymptotic 
normalization ratio f] = Ad/ As- Therefore the contributions from explicit two-body currents 
or from isobar configurations in the deuteron wave function should be relatively small at 

Encouraged by this result, in the present paper we study, within the framework of 
FFHD and using our Poincare covariant current operator, the deuteron form factors at 
7^ and in particular the effects produced by : i) different N — N interactions, and 
ii) different nucleon form factors models. We will also investigate the possibility to gain 
information from elastic e — d scattering on the neutron em structure, and in particular on 
the neutron charge form factor. Our preliminary results were already published in Ref. 

The plan of the paper is the following: in Sect. 2 the definition of our covariant 
current operator is recalled; in Sect. 3 the elastic deuteron form factors are expressed in 
terms of the matrix elements of the free current in the Breit frame; in Sect. 4 the front-form 
deuteron wave function and the explicit expressions of the current matrix elements in terms 
of the deuteron wave function are presented; in Sect. 5 our results on the dependence of 
deuteron form factors upon N — N interactions and nucleon em form factors are discussed 
and, eventually, in Sect. 6 our conclusions are drawn. 
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2 A covariant current operator within the front-form 
dynamics 



In this section we give the essential hnes for the definition of a current which satisfies ex- 
tended Poincare covariance, Hermiticity, current conservation, and charge normahzation, to 
be apphed to the calculation of elastic em form factors. 

Let us first consider the extended Poincare covariance. If the current operator J^{x) 
is defined in terms of J^(0) 

Jt'(x) = exp{iPx)J''{0)exp{-iPx), (1) 
then the Poincare covariance of J^{x) takes place if 

u(ir\j^m{i) = Lmj^io), (2) 

where L{1) is the element of the Lorentz group corresponding to / G SL{2,C) and U{1) is 
the unitary representation operator corresponding to / (see, e.g., [0). 

For systems of interacting particles the operator U{1) in general does depend on the 
interaction, and it is not trivial to build up a current which satisfies Eq. (^). Indeed, in 
order to fulfill this requirement the current operator has to be interaction dependent. The 
key property of our procedure |]10| for the definition of a Poincare covariant current operator 



is the following spectral decomposition of the current: 



j^(o) = ^n,j^(o)n,- (3) 

In Eq. (|^) Ilj is the orthogonal projector onto the subspace Hi = HiH corresponding to the 
mass Mi, the spin Si, and a definite parity, with H being the space of states describing the 
interacting particle system. This decomposition allows one to express the possible current 
operator dependence on the interaction as a dependence on mass and spin of the interacting 
particle system. 

In the FFHD, the seven Poincare generators belonging to the subgroup which leaves 
invariant the hyperplane x'^ = are kinematical. Then, as already mentioned in the intro- 
duction, the state of a system, \P,x)^ factorizes in a total momentum eigenstate, |Pi,P+), 
times an intrinsic eigenstate, \x) '■ 

\P,X) = \P±,P-')\X) , (4) 

In Eq. (i P+ = (Po + Pz)/V2 = p+ + ... + p+ and P^ = (P,., Py) =p^^ + ...+ p^i. are the 
plus and _L components of the total momentum, with pi, ...,Pn the individual momenta of 
the particles in the system. Because of the decomposition of Eq. (|^), the operator J'^(O) is 
fully defined by the set of matrix elements between initial, \P±,Pj'), and final, Pj"*"), 
total momentum eigenstates 

np:;Pj) = (Pl,i^'+|n,j^(o)n,|p^,p,+) (5) 
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The matrix elements between total momentum eigenstates, J^{Pl', Pj), correspond to definite 
values of masses, spins and parity, and are operators in the space Hint of intrinsic states. 
Through proper unitary transformations, the current operator J'^{Pl,Pj) in any reference 
frame can be defined in terms of the auxiliary current operators 

f{Ke,;M„M^) = {K[^ = 0, K'r^\U,J^^{0)U,\K,^ = 0,K+) (6) 

in the special Breit frame where the total three-momenta of the system in the initial state, 
Kj = —Kszi cind in the final state, K[ = Ke^, are directed along the spin quantization 
axis, z. In Eq. the initial and final plus components of the total momentum are 

Kl = ^[{M, + K'y/' - K], K'^ = ^[{M^ + K'y/' + K], (7) 

while K = Q/2, and q = K[ — Kj. It has been shown ||10| that the operator J'^(O) fulfills Eq. 
(0), i.e. is Lorentz covariant, if the current operators j'^{Kez', Mi, Mj) in the above special 
Breit frame are covariant with respect to rotations around the z axis. 

Since in the front form the rotations around the z axis are interaction free, the con- 
tinuous Lorentz transformations constrain the current j'^{Kez; Mi, Mj) for an interacting 
system in the same way as in the non-interacting case. The same property holds for the 
covariance with respect to a reflection of the y axis, Vy, and with respect to the product of 
parity and time reversal, 6, which leave the light cone = invariant, and therefore are 
kinematical. The full space reflection is the product of Vy and a dynamical rotation around 
the y axis by vr, while T = 6V, and therefore parity and time reversal do not contain an 
interaction dependence different from the one implied by rotations around y axis. As a con- 
sequence, the current operator satisfies V and T covariance, if it satisfies Poincare covariance 



and covariance with respect to Py and 9 ||10 



In conclusion, since in our Breit frame the extended Poincare covariance constraints 
for the auxiliary operators are the same for a non-interacting and an interacting system, the 
extended Poincare covariance is satisfied for an interacting system by a current composed 
by the sum of free, one-body currents, viz. 

j;_(ire,;M„M,) ^ (0, ir;+|n, j;_(o)n,|o, iT/) (s) 

where Jfred^) = Y.iLijfree,v with N the number of constituents in the system. 

In the elastic case considered in this paper (Mj = Mj = M; Si = Sj = S), the 
property of Hermiticity for the auxiliary operators reads as follows 

f{-K- M, M) = r{K; M, M)* , (9) 

where * means the Hermitian conjugation in the internal space Tiint- For =0 the 
property of Hermiticity reads j'^(0; M, M) = j^{0] M, M)* , while for | 7^ it becomes a 
non-trivial constraint and is satisfied if 

f{Kez-M,My = 

L[r,{-Ti)Y^D'[r^{-'n)]f{Kez- M, M)D'[r,{-'n)]r\ (10) 
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where D'^{u) is the matrix of the unitary irreducible representation of the group SU{2) with 
spin s, corresponding to m G SU{2), and r^l—w) represents a rotation by — vr around the x 
axis, i.e. D^[rx{—7i)] = exp^mSx) [|lOi . 

Let n be the projector onto the subspace of bound states |x) of mass M and spin 
S, and let J'^{Kez] M, M) be a current which fulfills extended Poincare covariance. Then 
a choice for the current compatible with the Hermiticity condition, Eq. (Il^), and with the 



extended Poincare covariance is 10 



r{Ke,; M, M) = ]^{J^{Ke,- M, M) + 
L^:[r,.(-7r)] exj9(m^,.) \r{Ke,-M,M)X exp(-mS,)}. (11) 

The second term in Eq. (|TT|), which ensures Hermiticity, introduces implicitely two-body 
terms in the current, because of the presence of the x component of the front-form spin 
operator, Sx- 

This current fulfills also the current conservation, which in the elastic case reads 

j^(ire,;M,M) = j+(ire,;M,M) (12) 
Indeed, as shown in ref. ITIl, in the elastic case the extended Poincare covariance and 



Hermiticity imply Eq. (|T2D, i.e., impose current conservation. 

In Eq. (|TT1) one has to choose a specific definition for the operator J^i^Kez] M). 
Unfortunately, one cannot simply adopt Eq. (^, because of the charge normalization con- 
dition, which implies 



j+(0; M, M) = ^iJ^iO; M) + J-{Q- M,M)} = v^eMH. (13) 

where e is the total electric charge of the system. Indeed, while the charge normaliza- 
tion condition is fulfilled by j^^^^{0; M, M), Eq. (|13D is not satisfied by |(j;^^,g(0; M, M) + 
jy^gg(0; M, M)). However, a possible choice is the following one: 

J+{Ke,;M,M) = (0, ir'+|HJ+,,(0)H|0, 
J^{Ke,;M,M) = (0, K'+|HJx/.ee(0)H|0, 7^+) 

J-{Ke,-M,M) = J+{Ke,-M,M). (14) 

The previous definition of the "— " component of JT^ is essential for the proper charge 
normalization of j^{Kez; M, M), because of the second term in Eqs. ( pT| , p!3D . 

In the elastic case, only 25" + 1 non-zero matrix elements of the em current defined 
by Eqs. (|ll],|14D are independent, corresponding to the 25* + 1 elastic form factors. The 
independent matrix elements can be chosen as the diagonal matrix elements of with 
S, >0 and the matrix elements {MSS^lj^iKe,; M, M)\MSSz - 1) of j^; with S, > +1/2 



T0| . Obviously, any other choice of the independent matrix elements to be used in the 
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calculation of the elastic form factors is completely equivalent, i.e. it will yield exactly the 
same results. One can immediately obtain that 



{MSS,\j+{Ke,;M,M)\MSS,) = {MSS,\J+{Ke,; M, M)\MSS,), (15) 
{MSS,\UKe,;M,M)\MSS:) = ^[{MSS,\J,{Ke,;M,M)\MSS',) - 
{MSS'JJ,iKe,;M,M)\MSS,)] (16) 

and therefore the elastic form factors can be evaluated in terms of the matrix elements of 
the free current only. It has to be noted that the matrix elements of both j"*" and jx have 
been shown to be real fiy . 



In the deuteron case, since 5" = 1, three matrix elements of the current are needed. 

3 Deuteron electromagnetic form factors 

The form factors A(Q^) and B{Q'^)^ which appear in the unpolarized cross section, and the 
tensor polarization, T2o(Q^), can be expressed in terms of the charge, GciQ"^), quadrupole, 
Gq^Q"^), and magnetic, Gm{Q^), elastic form factors : 



A{Q') = Gl + Ir^Gl + \rGl 



B{Q') = \r{l + r)Gl, 

T(n^,- ^V2[lrGl + AGQGc+2Gl] 

where r = (5^/(4m^), Q'^ = —qf^, rrid is the deuteron mass and / = 1/2 + (1 + r) tan^(6'/2), 
with the following normalization for the form factors: Gc'{0) = 1, Gq{0) = m^Qd, and 
Gm{0) = Hd'md/mp {nip is the proton mass). 

For the deuteron, the matrix elements of the current are related to the form factors 
GciQ"^), Gm{Q'^), Gq{Q'^) by the following general expression of the macroscopic current 
for spin 1 systems (as the deuteron) ||2^ 

2 

-{Gc - -^rGQ)gaf3 



+ Gm [g^qp - g^qc 



-Q^[Gc-{l + \r)GQ-GM]qaqp 

where \mdlSz) is the deuteron intrinsic eigenstate, ga(3 the metric tensor, es^ and e^, are 
the initial and final deuteron polarization vectors, respectively, (see Appendix A) and = 
\/2md\/l + r. 
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In FFHD, hadron form factors are often calculated in the reference frame where = 
0. If A and A' are the helicities in the initial and final states, respectively, and lyx = 
(A'l J"'"(0)|A), then, because of Hermiticity, V and T covariance, and covariance for rotations 
about the z axis, all the matrix elements lyx for the deuteron can be expressed in terms of 
hi, 



Iqo, Jio and h-i- As shown, e.g., in Refs. 
"angular condition", must be fulfilled in the q 



T^, H, the following constraint, usually called 
= frame, viz. 



:i + 2r)/n-/oo-(8r)i/2j^o + /i, 



0. 



(19) 



However, this constraint, which is related to the rotational covariance of the current, is not 
satisfied if the matrix elements lyx are calculated with the free operator, Jf'reei^) the 
g"*" = frame. Then, three out of the four matrix elements are usually defined through the 
free operator, while the fourth one is defined by Eq. (|19|). However, different choices of the 
three matrix elements to be calculated by the free operator are possible and therefore different 
prescriptions can be used to calculate the three physical form factors. As a consequence, 
within this approach there is a large ambiguity in the theoretical results (see 
furthermore 



20, 21 



e.g. 



18, 19 



3|), and, furthermore, different definitions are used for different matrix elements of 
the current. 

A relevant result of our approach is that, using in the left hand side of Eq. (18) the 



microscopic current defined by Eqs. (^,0), the extraction of elastic em form factors is no 
more plagued by the ambiguities, which are present when the free current is used in the 
reference frame where g"*" = 0. Indeed, using our current operator, it turns out that only 
three matrix elements jg, are independent, corresponding to the three elastic em form 
factors. For instance, one can consider the matrix elements Jq q, jii, jf g, which have been 
shown to be real ||T^. On the contrary, using the one-body current in the = frame, one 
has four independent matrix elements [|TH|. 

The form factors Gc, Gm, and Gq can be easily obtained from the matrix elements 
of the current in our Breit frame, since from Eq. (113) one has 



{mdll\j^{Ke;,;md,md)\mdll) = C ^ 
{mdlO\j^ {Ke^; rrid, md)\mdlO) = ('^ 
{mdll\j:^{Ke^;md,md)\mdlO) = C^r^'G 



Gc - 

Gc + ^tGc 



M 



(20) 



By means of Eq. 



2^ and using the properties ([TsD, (|l^) of the matrix elements 
jg, g^, the form factors Gc, Gm, and Gq can be expressed in terms of the matrix elements 
J't,s, = {mdlS^\J+{Ke^,md,md)\mdlS;,) and J|, 5^ = {mdlS'JJ^{Ke^,md,md)\mdlS^), 
i.e. in terms of the matrix elements of the free current, calculated in the Breit frame where 



One obtains 



the momentum transfer is along the spin quantization axis, z 

Gc = (2 J-+ + JoloKA G'm = (J?o-^o!i)C/(2v^), Gq = {J^^, - J,^,)C/{2r). (21) 
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Then, the deuteron magnetic moment, in nuclear magnetons, is given by 



771 1 

while the deuteron quadrupole moment is 

We stress that, as was shown in using the free current in the frame where = 0, in the 
hmit Q"^ ^ the angular condition is satisfied at the first order in Q, but it is violated at 
the second order. Therefore the angular condition is not a problem for the calculation of fid, 
while the quadrupole moment is not uniquely determined within that approach. 

From Eqs. (p^Tjpll) it is straightforward to obtain the expressions for the elastic struc- 
ture functions A{Q'^), B{Q'^) and for the tensor polarization T2o(Q^) terms of the matrix 
elements of the free current Js^^s^ ^s'^^s^'- 

Am = y [(Jolo)' + nJi^i)' + iJi,o - J.^f/'A 



6m 



^''^^ ^ - T A + B t^n\e/2) 



4 Deuteron front-form wave function and matrix ele- 
ments of the current operator 

We consider the deuteron as a system of two different, interacting particles with the same 
mass, m = (mp + m„)/2 (m„ is the neutron mass), and spin 1/2. For a system of N particles 
with four-momenta Pi {i = 1,2, N), FFHD internal variables ki, k^- can be defined, 
such that J2f=i ki = . The intrinsic three-momentum ki is the spatial part of the four- vector 

h = L[P{G)]-'pi, (25) 

where G = Pq/Mq is the four- velocity, and Pq = pi + ... + the total four-momentum 
of a system of free particles, with Mq = |Po| = I-Pq^I^^^- The matrix P{G) G SL(2,C) (see 
Appendix B) represents a front-form boost. The action of the boost L[(3{G)]~^ is such that 
Po = L[l3{G)]-'Po = [P'o^ = 0, Po+ = Mo, Pi = Mo]. 

Then the wave function for the deuteron internal state \mdlSz) = can be 



written as follows [30 



Xi,sAk±i^i(^u(^2) = {k,ai,a2\xi,sj = {k,(Ti,a2\R ^\^d)uj{ky^^, (26) 
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where ^ = pf /P~^, and k± = pi± — ^P±. The internal three-momentum is k = {k±,kzj, 
where k^ = (2^ — l)to{k), uj{k) = {m? + P)-^/^, and k = \k\. It can be easily shown that 
Mq = 2uj{k). The normalization of {k,(Ji,a2\xi,S:,) is such that 

E/l(*>..-.lx...)P(^=l (27) 
The matrix R is given by 

R = v{k,si)v{-k,S2) (28) 
where w(/c, s) is the Melosh matrix ||3l], |15| 

^ _ i^{k) +m + kz + ij^xky - dyk^) 

[2{uj{k)+m){uj{k) + k,)Y'^ ' ^ ^ 

while si, and S2 are the usual nucleon spin operators, ai and a2 the eigenvalues of si^ and 
S2zi respectively, and ai the Pauli matrix operators. The generalized Melosh matrix can also 
be written as 

v{—k,s) = exp{-ipn&), (30) 
with n = a fc)/A;_|_, by defining the angle (p 

(f = 2arctan — — , (31) 

Lj[k) + m — kz 

The angle (p will be used in the Appendix for the calculation of the deuteron quadrupole 
moment. 

The wave function for the deuteron internal state obeys the mass equation 

M^XiA ik±, ^, cTi, 0-2) = mlxi,s, ik±, ^, ai, (X2) (32) 



while the wave function in Eq. ( P6D is the usual solution of the "nonrelativistic" 
Schroedinger equation. Indeed, if in the front-form dynamics the mass operator M for the 
function is defined by = RM^R~^ = Mq + V with V the interaction operator, then 
the mass equation M^'^a = ''^l^d has the same form as the "nonrelativistic" Schroedinger 



equation in momentum representation 32, 15 



i.2 ^ 

(— + V)^rf(fc, di, as) = ^d^d(fc, ai, as) (33) 
m 

where 

V = V/Am, Ed = {ml - An?) /Am = + (4m) (34) 

with md = 2m + e^. Therefore the eigenvalue Ed of Eq. ( p3D can be identified with the 
deuteron energy e^, if the small quantity ed/{Am) is disregarded. It has to be noted that, 
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in the case of the N — N interactions of the Nijmegen group Ed is directly linked 
through Eq. ( P^D to the deuteron energy used in their fits. For the continuous part of the 
two-nucleon spectrum the mass equation is identical to the " nonrelativistic" Schroedinger 
equation in momentum representation |]I6| . Therefore the operator V has to satisfy the same 
constraints of the potential as in nonrelativistic quantum mechanics and can be chosen to 
have any of the forms usually employed for the N — N interaction in nonrelativistic nuclear 
physics. 

Since the wave function is an eigenstate of the standard nonrelativistic spin oper- 
ator i, |16l 

Snr = l{k) + S1 + S2 (35) 

where l{k) is the usual orbital angular momentum, the Clebsh-Gordan coupling coefficients 
can be used. Then the internal deuteron wave function Xi,Sz(^? cti^ ^"2) with polarization 
vector 65^ (see Appendix A), is given by (cf. W) 



1 3kk- 
Mk)Sij - -^i^ij p^)v'2(A:) 



2^7rMk){esJ, 



Xoik)5^j + ^J^)Mk) 



v{k,si) ^aiayv{-k,S2)* 



0-1,0-2 



(36) 



where a sum over the repeated indices i,j = 1,2,3 is assumed and Xo{k) = ipo{k) — 
{1/V2)^2{k). The wave functions fo{k) and (p2{k) coincide with the nonrelativistic S and 
D waves in momentum representation ||32[. The normalization of fo{k) and ^2{k) is such 
that /[(y9o(^)^ + ^2{ky']dk = 1. For the calculation of the matrix elements of the current it 
will be useful to put the internal deuteron wavefunction in a more compact form 



Xi,sAk,a,,a2) = 2^7TMk){Ss.)jF,,{k) [a{k) + za ■ A(fc)],^,„ 



(37) 



where 



Xo{k)6ij + -^y^Mk) 



(3J 



and 



Ci{k) + la- ■ Di{k) = v{k, si) ^aiay[v{-k, 32)]" 



1,2,3. 



(39) 



In this paper the matrices Ci{k) + t& ■ Di{k) will be called "generalized Melosh matrices for 
the deuteron wave function". Explicit expressions for the real quantities Ci{k), Di{k) can 
be found in Appendix C. 
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The matrix elements J^s' ^e easily calculated, by using the action of the free 
current on a two-body state \P±, P~^)\xs,sJ pOj : 



w{pi,ai){k, (Ti,a2\xs,sJ 



f:{{p'i-Pif)-mp\~Pif) 
1 



{pi+PiY 



(40) 



where, in our case. 



^/^ee(0) = ^;(0) + JaO). (41) 

In Eq. (^OD w{p, a) is the front-form Dirac spinor (see Appendix B), while = f^+fe 
and /m = /m + fm the isoscalar electric and magnetic Sachs form factors of the nucleon. 

An explicit calculation, with the help of the matrix elements of the 7 matrices between 
front-form Dirac spinors reported in Appendix B, shows that, as a consequence of Eqs. 



V2md Yl J Xi,s.(^', 0-^,0-2) 

dk' 



a?rn? + b'^k'^. 



f 



Xl,5.(^,C^l,0^2)(^C 



Al/2_ 



{27ifuj{k% 



(42) 



(Xi,sdX(i^e-;,m,,m,)|xi,s.) = {xi,s'MK'+\rf^^M\^,K+)\xi,s^ 

\AmkM' - r^)[am + ih{ak)^] 



Y: J Xi,s'Sk',a[,a,rl 

Xi,5,(A;, ai,(T2 
where 



+ fmV^^x + ih{may + kyd^)] 



dk' 

(2^0M^ 



(43) 





1/2 






1/2 


h = 




1/2 






+ 

















nl/2 



(44) 



and the form factors /^^ and are functions of {p'l — piY- In our Breit reference frame, 
where i^j, = and gl = 0, the relations between the internal (/c_l, kz) and individual nucleon 
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— * — * 

variables, in the initial, xi ^^(fc, ai, (T2), and final, xi s' {k' , cr[, , wave functions are given 

by 



/+ 

Pi 

K'+ 



1 + - = / . = -(45) 



K' 



m 



with K = K/irtfi. It is important to note that nucleon form factors cannot be factorized out 
in the current matrix elements, since from Eq. (^3p one has 

(K - Pif = -^r{m' + kl)/m ^ -Q'. (46) 
By using the expression (pUf) of the deuteron wave function, a direct calculation shows that 



{XhS,\J^{Ke^,md,md)\xi,s.) = V2md{esjy{esjj 



u;{k%^ 



nl/2 



XQ{k')5i>j, + '-^^V2{k') 



Xoik)5,j + ^J^Mk) 



-Tr I [v{-k', S2)faycri,v{k', si) 
v{k, Siy^cricry[v{-k, S2)]* \ dk' 



am{f'/ - f^) [am + tb{d-k)±] 
a?rn? + 



(47) 



{Xi,s',\Jx{Ke^,md,md)\xi,s,) = ies'J*,ies,)^ 

/ 1 \ r- ^k-jk'j 



Xo{k)6i'j' + -^^V2[k) 



^Tr i [v{-k', S2)fayai>v{k', si) 



V2k 



I 


'uj{k)' 




_U}{k')_ 


\v2{k) 




- fin) [am 







1/2 



+ 



1 1 dk' 

f^[ak^ + ibimdy + kya^)\^v{k, si)~^aiay[v{-k, S2)]* j — 



(4J 



where the superscript T on a matrix indicates the transposition of the matrix and a sum over 
the repeated indices is understood. By means of the matrices Ci{k) + 1& ■ Di{k), 

equations (^7|), ( ^8l) can be rewritten as follows : 



{Xi,sJJ^iKe^,md,Tnd)\xi,sJ = y^rnd{esSj'i^sJ j J 
-Tr{[Ci>{k') - to- ■ D,,{k')] [A+ + la ■ B+] [Ci{k) + la ■ Di{k)]} dk' 



iAkW 



F,,,{k')Fi^{k) 



(49) 
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{Xi,si\JxiKe^,md,md)\xi,sJ = (e5^)*,(e5j 



/ 


'u;{k)' 




_uj{k')_ 



nl/2 



where 



1 fik' 
-Tr{[C,{k') - i& ■ D,{k')] [A, + i& ■ B,] [Q{k) + la ■ A(fc)]} — 
^ 4 



A + l(T ■ rS = — ttttt; h 



and 



Amk^iP/ - f^)[a'm + ib{ak)±] 



+ fmi'^^x + ih{may + kyO^)]. 



(50) 



(51) 



(52) 



It is straightforward to see that A^ is proportional to the quantity a, while B^ and B^. are 
proportional to 6, defined in Eq. (0). All the quantities A'^ ,B^ ,Ax,Bx.,Ci{k)^Di{k) are 
real. 

By an explicit calculation of the traces in Eqs. (^), (0) one has 



r \uo{k)f' 

{Xi,s,\J^{Ke^,md,md)\xi,s.) = ^/^md^esSfi^s-^j J '^g^ 

{A+ [Ce{k')Q{k) + Deik') ■ D,{k)] - 

B+ ■ [Ci,{k')Di{k) - Di,{k')Ci{k) - Di,{k') A A(fc)]} dk' 



1/2 



F,,,{k')F,^{k) 



(53) 



/ 


'u}{ky 




_U}{k')_ 



1/2 



{Xi,s'JJx{Ke;,,md,Tnd)\xi,s^ = ies'J*,ies, 
{A, [Qik')a{k) + A'(fc') • A(fc)] - 

rik' 

B, ■ [C,{k')Di{k) - D,ik')Q{k) - D,{k') A A(fc)]} — 



F,,,{k')F,,{k) ■ 



(54) 



It has to be noted that the integrals in Eqs. (p3| , 54) are real. Therefore, since the matrix 
elements Js^^s^ ^s'^^s^, (s^e the end of Sect. 2), only the real part of (65/ )*,(e5J^■ 

can contribute to these matrix elements. 



5 Numerical results for the deuteron form factors 
5.1 Deuteron magnetic and quadrupole moments 

The direct evaluation of magnetic and quadrupole moments through the limits of Eqs. (p^,p^) 
implies very delicate numerical problems and then a careful analitical reduction of these 
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equation is needed. For the sake of completeness we report in Appendix D the exphcit 
expressions that have actually been used. Magnetic and quadrupole moments have already 
been calculated in Ref. [0 for a variety of N — N interactions. In this paper we recall 
our main results, which are summarized in Table I. In the table the values of the magnetic 
and quadrupole moments calculated with many N — N interactions, already shown in Ref. 
||13| , are reported together with the values obtained using the local Nijmegen2 interaction. 



which was not considered in Ref. 13 



The standard non-relativistic results obtained with a one-body current crucially de- 
pend on the asymptotic normalization ratio t] = Ad/ As of D and S wave functions and 
on the D— state percentage in the deuteron, Pd, but one cannot obtain at the same time 
the experimental values for both and Qd- Using the free current within the FFHD in 
the g"*" = reference frame, the relativistic correction (RC) turned out to be very small for 
Qdi while for /i^ it could explain only part of the disagreement with the experimental value 
Q. On the contrary, in our Poincare covariant calculation the RC's bring both ^d and 



Qd closer to the experimental values, except for the charge-dependent Bonn interaction 
We wish to stress that our current operator and the one used in Ref. are different, since 
both of them are obtained from the free one, but in different reference frames, related by an 
interaction dependent rotation. As was already observed for the nonrelativistic calculations 
of Qd [il, |39|, we have shown in Ref. that a remarkable linear behaviour against the 



asymptotic normalization ratio, rj, holds for both the deuteron moments calculated within 
our approach (the values of Hd and Qd corresponding to the Nijmegen2 interaction obey 
precisely the same trend as the other interactions). The values of fj,d and Qd, suggested by 
this linear behaviour in correspondence of the experimental value of r] [rj^^P = 0.0256(4) ||40[| ) 
differ from the experimental ones (/i^ = 0.857406(1) [|T[ and Qd = 0.2859(3) g) only by 
0.5% and 2%, respectively, i.e. much less than for the non-relativistic results. The RC to 
lid is rather large and the total result becomes slightly greater than /i^^^, while the nonrel- 
ativistic one is smaller. This shows that, within our framework, even the sign of explicit 
contributions of two-body currents is different from the one needed in the non-relativistic 
case. In conclusion, it appears that, within our approach, the total contribution of two-body 
currents (from meson-exchange, Z-graphs, etc.) and isobar configurations has to be relatively 
small at Q'^ = 0. 



5.2 Deuteron form factors and N — N interactions 



Let us first compare in Figs. 1 and 2 our relativistic results for A^Q"^), B{Q'^) and T2o(Q^), 
obtained using the RSC interaction and the Gari-Kriimpelmann nucleon form factors 
with the corresponding nonrelativistic results. Following Lomon [24|, the latter ones 



42 



have been obtained by using the exact relativistic relations between the deuteron form factors 
and the current matrix elements, within the Breit reference frame where the momentum 
transfer is directed along the z axis but with nonrelativistic expressions for the matrix 



elements evaluated in impulse approximation p4 |. 

In order to have a closer insight to the form factor behaviour, in addition to the 
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usual plots for A(Q^) and B{Q'^) in a logarithmic scale, shown in Fig. 1, we report in 
Fig. 2(a) the quantity A^Q"^) divided by the factor (Gf, ■ F), with = (1 + gVO.Tl)-^ 
and F = (1 + Q^/O.l)"^-^, in a linear scale, and in Fig. 2(b) the quantity Tm{Q'^) = 
[GM{Q^)mp/{fXdmd)]'^ divided by the factor (G^-Fi), with Fi = (l + QVO.l)-^. As it is clear 
from Figs. 1 and 2, the differences between relativistic and nonrelativistic results are a few 
percent for < 0.1[GeV/c)'^, while become large as increases. For A(Q^) the differences 
are larger than 20% already at > 0.2{GeV/c)'^ and are of orders of magnitude for > 
2{GeV / cY . For B{Q'^) the relativistic and nonrelativistic results differ by 50 — 100% for 
> 0.3(Gel^/c)2, while for T2o{Q'^) they considerably differ for > 0.5(^61^/0)2. Yigs. 
1, 2 we have also reported by dashed lines the results obtained by keeping fixed the argument 
of the nucleon form factors in Eqs. ( ^ , ^8] ). The effects of factorization become large for 
/l(Q2) and S(Q2) at > i[GeV/cf, while for T2o{Q'^) already at > Q.^[GeV/cf. From 
Fig. 1 it appears that the nonrelativistic approach is able to give an overall description of 
the data for A^Q"^), B{Q^), and T2o{Q'^). However, this description is not accurate, even at 
very low values of the momentum transfer, as one can see in Fig. 2 (a) and, furthermore, it 
strongly depends on the N — N interaction and the nucleon form factor model. For instance, 
using the CD — Bonn interaction |37] and the nucleon form factors by Hoehler et al. [ |60[] , 
for v4(Q2) and T2o{Q^) the agreement is completely lost at > OAlGeV/c)"^. 

A comparison of our results with the deuteron form factors obtained by using the 
same N — N interactions and the same nucleon form factors, but within different relativistic 
approaches, for instance within the front-form calculation of Ref . , can also be interesting. 
Using the Paris interaction and the form factors of Ref. |Q, large differences have 
been found for A^Q"^) at > 2{GeV/cy, which become of orders of magnitude at = 
Q{GeV/cf (see Ref. |27|]). For B{Q^) we found a minimum around = 1.8{GeV/c)'^ 
instead of = l.6{GeV/c)'^ as in Ref. 0, and for T2oiQ^) a zero at Q"^ = lAiGeVjcf 
instead of \.2{GeVlcf. 

The results obtained within our approach with different N — N interactions are an- 
alyzed in Figs. 3 and 4, using the nucleon form factor model by Hoehler et al. [Q. We 
consider the old RSC interaction [^| and recent realistic interactions, able to describe the 
two-body data with a reduced ~ 1- Iii particular we study the A^IS interaction by the 



Argonne group [^, some interactions by the Nijmegen group {Nijmegenl, Nijmegen2, 
Nijmegen93, Reid93) [^, and the charge- dependent CD — Bonn interaction by the Bonn 
group ||3^. The results for the Reid93 interaction are essentially equal to the results of the 
^41/ 18 interaction and are not reported in the figures. 

The effects of different interactions are large for ^(Q^) at > l{GeV/cy, while for 
B{Q^) and T2oiQ'^) already at > 0.5{GeV/cf. It can be noted that the CD - Bonn 
interaction, which is characterized by a larger non-locality, yields larger differences with 
respect to the other interactions. At low values of (Q^ < OA{GeV/c)'^), where the nucleon 
form factors are better known, a simultaneous description of the experimental data for AIQ"^), 
B[Q'^) and T2q{Q'^) is achieved. The dependence on the nucleon-nucleon interaction in this 
region is minor, although not negligible (see, in particular. Fig. 4 (b)). 

For the mentioned interactions and using the Gari-Kriimpelmann nucleon form factors 
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42| , we report in Fig. 5 (a) the value of corresponding to the minimum of B{Q'^) and 
in Fig. 5 (b) the value of corresponding to the second zero of T2q{Q'^) against the 
nonrelativistic S-state kinetic energy, T^, in order to find a correlation between different 
effects of the N — N interactions. For both quantities a distinct linear behaviour is clear: a 
lower value of Ts yields a minimum for B{Q'^) and a zero for T2o(Q^) at a larger momentum 
transfer. Analogous results can be obtained with different nucleon form factors, as the ones 
of Ref. From Figs. 3(a), 4(a) it is clear that for > liGeV/cY' a similar correlation 
holds for A^Q"^)^ i.e. a lower value of T5 yields a lower value of A{Q'^). It has also to be 
noted that the AV^18 and Reid93 interactions, which give essentially the same results for 
Ai^Q"^), BiQ"^) and r2o(Q^), have the same S-state kinetic energy. 

Let us note that recent measurements of the S — D mixing parameter, ei, point to 
a stronger tensor force than the one exhibited by the interaction models we have analyzed 

pT| . In turn, a stronger tensor force is favoured by a high degree of locality, which yields 



significantly larger kinetic energies and, in particular, larger values of T5 Then, by an 



extrapolation of the linear relations found above, one can argue that a N — N interaction able 
to reproduce these recent measurements of ei could yield, on one hand, agreement between 
experimental and theoretical values for T2o(Q^) and, on the other one, a minimum for -B(Q^) 
slightly lower than the value indicated by the available experimental data (around = 
1.6{GeV/c)'^ instead of = LSlGeV/c)"^). Therefore, if new, more precise experimental 
data for B{Q'^) will show such a lower value for the position of the minimum, both BiQ"^) 
and (Q^) could be reproduced by a novel N — N interaction, without a relevant role for 
explicit two-body currents. 

5.3 Deuteron form factors and nucleon electromagnetic form fac- 
tors 

In order to investigate the effects of the nucleon form factors on the deuteron form factors, 
we have displayed in Fig. 6 our results obtained with the Nijmegen2 nucleon-nucleon 
interaction and corresponding to the nucleon form factor models of Refs. ^2|], |Q, and 



For A(Q^) the differences between different models are very large at > 0.5(GeV/c)^, 
increase as increases, and can be related to the sizeably different behaviour of G^^Q"^) 
for the various models. The influence of the nucleon form factor models is less marked in 
B{Q'^), while, as already known [Q, the tensor polarization is essentially independent of the 
nucleon form factors. 

Therefore, the linear behaviour of the locations of the minimum of B{Q^) and the 
second zero of T2o(Q^) vs. T5 is substantially independent of the form factor models, as 
well as the conjecture at the end of the previous paragraph. As far as AIQ"^) is concerned, 
one could try to exploit the strong dependence of A{Q'^) on G^{Q'^) to gain information on 
G^iQ"^) by a fit of the A^Q"^) experimental data, following a procedure analogous to the one 
used, in a nonrelativistic context, by Platchkov et al. [0. Obviously the results of this fit 
will be different for different interactions. Another possibility to be studied in our covariant 
framework is obviously the role of isobar configurations in the deuteron state (see, e.g.0) 
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and of explicit two-body contributions in the e.m. current (see, e.g. 0). As already noted 
T(]|,[Bn], these contributions have to be Poincare covariant, and to satisfy Hermiticity and 



current conservation by themselves. We intend to perform such a fit and to study these 
contributions elsewhere. 



6 Conclusions 

In this paper the deuteron form factors A(Q^) and S(Q^), and the tensor polarization T2o(Q^) 
have been evaluated in the framework of front-form Hamiltonian dynamics, using a Poincare 
covariant current operator, without any ambiguity. The current is built up from the free one 
in the Breit reference frame where q is along the z axis and fulfills parity and time reversal 
covariance, as well as Hermiticity and current conservation. 

Large differences have been found between the results of calculations performed within 
a nonrelativistic framework and within our Poincare covariant approach. These differences 
become huge at high momentum transfer, as expected, but are relevant for accurate calcula- 
tions even in the limit of zero momentum transfer, as is clear from our results for the deuteron 



magnetic and quadrupole moments |]I3| . Large differences have also been found with respect 
to a front-form approach which ensures Poincare covariance by different definitions for differ- 
ent matrix elements of the current operator (0). Our current operator, which was already 
shown to be able to describe the deuteron magnetic and quadrupole moments, is also able to 
simultaneously reproduce the three deuteron form factors at low momentum transfer, where 
the nucleon form factor are better known and the effects of different interactions are minor. 

The effects on the deuteron form factors of different nucleon-nucleon interactions and 
different nucleon form factor models have been studied. The different nucleon form factor 
models strongly affect y4(Q^), while the different interactions have large effects on A,B 
and T20. These effects are linked to the S'-state kinetic energy in the deuteron, which, in 
turn, is related to the degree of non-locality of the interactions and to the strenght of the 
tensor force. A novel N — N interaction with a strong tensor force, able to reproduce the 
recent measurements of ei, would be helpful to describe the deuteron form factors and, in 
particular, to offer a solid ground for the study of the neutron charge form factor from the 
analysis of A{Q'^). We stress the relevance of a well defined relativistic approach to gain 
reliable information on the nucleon-nucleon interaction and the nucleon form factors. 
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In this Appendix we report for the sake of completeness the expUcit expressions of some 
useful quantities needed for the calculation of the deuteron em form factors. 



A Polarization vectors 

The deuteron polarization four- vectors, e^^, in any reference frame can be obtained by a 
proper boost from the polarization vectors in the deuteron rest frame, es^{rf) = (e^j = 
0,65 J, with 



e+i = (1, I, 0), e_i = 4^ (1, -I, 0), eo = (0, 0, 1). 



V2 



V2 



(55) 



In our Breit frame, where P± — = 0, the transverse deuteron polarization vectors, in 
both the initial and final states, read as follows: 



e±i = e±i = (0, 1, ±t, 0), 

while the longitudinal polarization vector in the initial state is 

1 



and in the final state is 



eo = —{-K, 0, 0, ^Jmj + K^) 
ma 



e; = ^(i^, 0, 0, ^Jml + K^). 



(56) 



(57) 



(58) 



B Front-form Dirac spinors and matrix elements of 7 
matrices 



Adopting the following representation for the 7 matrices 

7° 






1 




1 




-ai 


1 







-1 


Oi 



(59) 



with i = 1, 2, 3 and the Pauli matrices, the front-form Dirac spinor ^(p, cr) can be written 
as 



w{p^ cr) = \fm 



(60) 



where x((7) is the ordinary spin 1/2 spinor describing the state with spin projection on the 
z axis equal to a and the matrix ^{g) has the components 



(61) 
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with g = p/m. 

One can immediately obtain 

w{p', (r')w{p, a) = -^=^==(a'| [m(p+ + - ta^^p-^Py - p'^Py) + idy{p-^p'^ - p'+p^)] |(t), 

(62) 

with normahzation 

(63) 



w{p,a')w{p,a) = —-{(j'\m2p'^\(T) = 2m5„(ji. 

p+ 

The matrix elements of the 7 matrices, needed for the calculation of the deuteron form 
factors, are 



w{p' , a')^'^w{p, a) = 2a/p+p'+5o 



(64) 



w{p', a')-f^w{p, a) = , ] {a'\ imq+ay + p^p'^ + p'^p^ + ia^{p'+py - p+p' ) |cr). (65) 
In our special Breit frame Eqs. (^) and (|65D become : 

=(cr'| m(p+ + p'+) + iq^{cr^ky - aj^fc^.) |cr) 



w(p',a')w(p,o-) 
w{p',a')'j^w{p,a) = 



= ((t'| imq^dy + (p+ +p'+)A;x. + lo^kyq^ \a). 



(66) 
(67) 



C Generalized Melosh matrices for the deuteron wave 
function 



The generalized Melosh matrices for the deuteron wave function have been defined in Sect 
4 as the matrices 

Ci{k) + id- ■ Di{k) = v{k, Siy^aiay[v{-k, S2)]*, i = 1,2, 3. 

From the expression (^) for the matrix v{k, s) one obtains 



a{k)=Af 



S2i'm + 



63,171 + 



[D,{k%=N 

[D,{k)]y=Ar{e,Ak) 
[D,{k%=N 



m + uj{k) 

zk'i — osi 



kyki S^ik 



m + Lu{k) 



Oii-m + 



where 



m + uj{k) 



1 



MoJai-0 Jm^ + k 



(6J 



(69) 

(70) 
(71) 
(72) 
(73) 

(74) 
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D Deuteron magnetic and quadrupole moments 



In this Appendix we illustrate the main steps for the calculations of the deuteron magnetic 
and quadrupole moments from Eqs. (^,^). To this end, expansions in k = i/r = Q/{2md) 
of the quantities a and b (Eq. (^4|)) up to the first order 

a = 2, b = 2^ (75) 
and of the quantities ^ and A;^ (Eq. (^S])) up to the second order 



e = r - 2k{1 - C) - 2«^(1 - a h = K- + ^(^')|^(4e' - 3) (76) 

will be needed, since the intrinsic moment in the final state, k', is the integration variable in 
the integrals for the calculation of the current matrix elements. 



D.l Magnetic moment 

The deuteron magnetic moment is given by Eq. (p^) 



(77) 



where the matrix elements J'^q and J'qi can be obtained by Eq. 



Let us preliminarly 



note that J^iq and J^q^ have the same expression, but for the exchange of the role of initial 
and final variables in Fi j{k), Fi/j/{k'), and in the quantity between curly brackets in Eq. 
([5^ ) (we recall that only the real part of (e+i)^ gives a non vanishing contribution to the 
matrix elements). In order to obtain the magnetic moment, one can expand [J^iq — J^oi] ^ ^ 
function of k, and consider only the terms which are linear in k (indeed, because of Eqs. (^) 
and (^61) , the current matrix elements are functions of k). As a first step, by using Eq. ([75|) 
and Eq. ([76|), we expand the quantity between curly brackets in Eq. ( p^ at the first order 
in K. We obtain a term independent of k and a term linear in k, which is identical, but with 
opposite signs, for the two matrix elements JT^q and J^qi- It is clear that in correspondence 
to the latter term one can evaluate the radial wave functions in Eq. (|5^ with the same 
argument k. After an integration over the polar angle [fc = {k, 6, 0)] one has: 



/.. = -lim "^^[^--^^] + 



Tcmrn 



- d{kl) ^f-Xo{k) 

Jo J-oo 



mpTT 
2mm3 



d{k 



''kl 



°° dk^ 



Xoik) 



Xo{k) 1 + 



kl 



2m{uj{k) + m) 



Xoik) + -^Mk) 



ris r 

J m J t 



3(p2{k) 

uj{k) 
uj{k) + m 



k^ + kl 
2V2k^ 



+ 



(78) 
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where the first term and the last two terms correspond to the zero and first order terms in 
the expansion of the curly bracket of Eq. (0), respectively. In Eq. ([75D, T is given by the 
following expression 



uj{k') 



^ = 37r/ d{kl)kll ^ 



nl/2 



V2k'' 



V2k^ 



3kz{k'^ + k^k'^ 



<P2{k)(f2ik') 



(79) 



and, according to the observation at the beginning of this subsection, JF' has the same 
expression, but for the exchange of k and k' in the quantity between square brackets. 
The limit in Eq. (^) can be easily handled and one obtains the final result 



mnip 

fid = Svr 2- 

mj JO 



k'^dk / dicosO) 
Jo 



[^[Mk)ni - cos'e) + f-xoik) 



k^ 

Xoik)^ 



Xo{k) + ^<^2(A;) 



J n 



Mk)r + ki] 

(m2 + kiy 

Xo{k) (2 

uj{k) 
uj{k) + m 



k' 



m{uj{k) + m) 



+ ^Mk) 



[1 + cos^e) 



V2 



+ 

BO) 



In Eqs. ( [78| , pOD the nucleon form factors and have to be evaluated in the limit 
Q ^ 0, i.e. /r(0) = 1, /-(O) = 0.8797. 

The nonrelativistic result for fi^ can be immediately recovered from Eq. ( ^0|) in the 
limit m 00. 



D.2 Quadrupole moment 

The quadrupole form factor (see Eq. (pTD) is given by 

V2md [JoTo - Ji:i] 



Gc 



vr+7 



(81) 



The proper combination (j7o\) — J7i\) of the matrix elements of J''^ can be directly calculated 
from Eq. ( |5B| ) by using Eqs. (|5BD, ( ^^ and the explicit expressions for the the quantities 
Ci{k), Di{k) given in Appendix C. One obtains 



Jo,o - J^A = mdV2 



f 


\^ik)^'] 


1/2 - 




Hk')^\ 





dk' 



(82) 



where 



E = -Xoik)xoik')[l-cosiif-ip')] + 
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+ 
+ 



and 



2FA;'2 ^2^' " + ^.^Dcos(<^ - + {k, - k',)k^sin{^ - ^')] (83) 



H = -Xo{k)xo{k')sin{ip - <^') + 



V2k'^ 
9ip2ik)Mk') 



1 3 

(kl - -k\)sin{(p - if') + -kzk±cos{(p - (f') + 

I 1 3 
(k'z - 2^1)sin{(f - If') - -k',k^cos{if - (f') + 



2k2k'^ ^^^^^ - 2^' ^f^^' + fc,A;:)s2n(<^ - <^') + {k, - k',)k^cos{^ - <^')] (84) 

The angle y?' is defined by Eq. (|3l|) with A; replaced by k' . 

The expression for Gq given by Eqs. ( pT| , p2| , p3| , p^ holds at any value of Q^. For the 
evaluation of the quadrupole moment 



lim — - 



/ 


W{k)i'] 


1/2 - 




[uj{k%\ 





"^-^^^^ a2m2 + 62A;2 " ^ 



dk' 



B5) 



an expansion of [j7o|o — Ji,i\ at the second order in k is needed. 
Let us note that at the first order in k one has 

, k^K 
^ ~^=e(a;(A;) + m) 

and, as a consequence, the quantity H is of the first order in k 

H = kHi + 0{k^) 

with 



(86) 



i7) 



2^ [u;(A;) + m 

' ^ ^^,2 ^^ 3v92(A:) 2V2yo(fc) + V^2(fc) + 3^2(A:) ^ ,V ' 



9uj{k) 
V2k^ 



k 



dk 



dk 
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Since b is also of the first order in k, (see Eq. ([75D), in Eq. (|85D one can take a = 2 and 
disregard b'^kj_ with respect to a^m^ in the hmit — » 0. As a result one has 



= Qdl + Qd2 







1/2 - 


f 







dk' 



^9) 



where 



Q 



dl 



Q 



d2 



— lim 



^0 g2 



'u{k)e 



1/2 



E dk' 



(90) 
(91) 



In the integral of Eq. (^) each quantity has been evaluated at = (i.e., ^ = k^ = k'^ 
and 



Eo = E{Q' = 0) = i-kl - kl 



3^2ik){2V2Mk) + V>2{k)) 



(92) 



To evaluate Qd we need an expansion of the integral in Eq. (pT|) up to the second 
order in Q. By using the expansions of ^ and k^ up to the second order in k given in Eq. 
(|76|), one obtains 



u{k%^ 



nl/2 



K 



1 + + —Vt2 



K 



where 



El 



4^-3 
2^' 

3cj(A;')A;; 

y/2^'k'^ 



E = Eq + kEi + ^-^2 



16^'2 - 36^' + 21 



(93) 
(94) 

(95) 



3P / 

'>^,(fc')U'o(^') + 



A; 



2v^ ; 



+ 



E'2 



kl 



dk' ' dk' ' 72 dk' 



(96) 



k'^^"^{uj{k')+my 



9cj(A;') 
V2k' 



uj{k') + m) 
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dk' 

3A;' 



k' 



+ 



dk' 



+ 



^ dk' 



(kj OMk') , ,j2 d'Mk') \ 



+ 



+ 



k± 
k' 



l^Mk')Mk'mJ - kl) - 6A::Vo(A:')^^ + ^[v^2(A:')]^+ 



Furthermore, because of Eq. i 

r:{{p\-Pif) = i + Q' 

where 

frln' = 6 



one has 



1 - 



J Q2=0 



"1 



2 I 2 

ep ' en 



JQ2=0 



(97) 



(9J 



(99) 



is the sum of the squares of the proton and neutron charge mean square radii (let us recall 
that {p'l —piY ^ 0). Then, since only the second order terms in the expansion of the integral 
in Eq. (^) can give a contribution to Q^, one obtains 



(12 



Ami 



Js\2 



dk 



(100) 



where each quantity has to be evaluated at = 0. 
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Interaction 



D 



1^ 



/id (LPS) 



Qd (LPS) fm' 



Exp 
RSC H 
Avl4 134 



Paris 

Avl8 12|] 
Nijm93 [H 
RSC93 

Nijml m 



Nijm2 |6 



CD-Bonn [|37 



6.47 
6.08 
5.77 
5.76 
5.75 
5.70 
5.66 
5.64 
4.83 



0.0256(4) [| 
0.0262 
0.0265 
0.0261 
0.0250 
0.0252 
0.0251 
0.0253 
0.0252 
0.0255 



0.8429 
0.8451 
0.8469 
0.8470 
0.8470 
0.8473 
0.8475 
0.8477 
0.8523 



0.857406(1) H 
0.8611 
0.8608 
0.8632 
0.8635 
0.8629 
0.8637 
0.8622 
0.8652 
0.8670 



0.2796 
0.2860 
0.2793 
0.2696 
0.2706 
0.2703 
0.2719 
0.2707 
0.2696 



0.2859(3) I? 
0.2852 
0.2907 
0.2841 
0.2744 
0.2750 
0.2750 
0.2758 
0.2756 
0.2729 



TABLE CAPTION 

Table I. Magnetic moment (in nuclear magnetons) and quadrupole moment for the 
deuteron, corresponding to different N — N interactions; /i^^ and Q^^ are the nonrelativistic 
results, fid (LPS) and Qd (LPS) our results; Pd is the D-state percentage, and r] = Ad/ As 



the symptotic normalization ratio (this table is taken from Ref. [|13|, a part from the results 
for the Nijmegen2 interaction). 



Table I P.M. LEV, E. PACE, G. SALME 



29 



FIGURE CAPTIONS 



FIG. 1. (a) Deuteron form factor v4((5^) obtained using the RSC N — N interac- 
tion and the Gari-Kriimpelmann nucleon form factors |4^. Sohd hne: full result of 
our approach with the Poincare covariant current operator. Dashed hne: the argument of 
the nucleon form factors, {p'^ — pi)^), is replaced by —Q^. Long-dashed line: nonrelativistic 
result obtained with exact relativistic relations between deuteron form factors and current 
matrix elements, within the Breit reference frame where q = |^ , but with nonrelativistic 
expressions for the matrix elements evaluated in impulse approximation |2^. Experimental 
data are from Ref. |Q (open squares), Ref. (triangles), Ref. |^ (diamonds), Ref. []TT| 
(full dots) and |jl2[ (open dots), (b) The same as in (a), but for B{Q'^). Experimental data 
are from Ref. (open dots), Ref. [^] (open squares), Ref. (full diamonds), Ref. pO | 
(triangles), Ref. ||5T| (full squares), and (open diamonds), (c) The same as in (a), but 
for r2o(Q^). Experimental data are from Ref. ||5^ (open dots), Ref. (full triangles), Ref. 
|]55| (open triangles), Ref. |Q (full dots), Ref. ||5^ (open squares), Ref. ||58[ (full squares), 
and Ref. [^] (diamonds). 



FIG. 2. (a) As in Fig. 1 (a), but for the reduced form factor A{Q'^)/{G']j ■ F) with 
Gd = {1 + gV0.71)-2 F = (1 + gV0.1)'2.5^ ^s in Fig. 1 (b), but for the reduced 
form factor rM(Q2)/(G'2^.Fi) with rM(g2) = [GM(g')mp/(/i,md)]2 and Fi = (l+gVO.l)-^. 
Experimental data are as in Fig. 1. 



FIG. 3. (a) The deuteron form factor AiQ"^) obtained using our Poincare covari- 
ant current operator, different N — N interactions and the nucleon form factors by Hohler 
al. 



et 
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Solid line: RSC interaction dashed line: Al^lS interaction p5[; dot- 

dashed line: Nijmegenl interaction; long-dashed line: Nijmegen2 interaction; short-dashed 
line: Nijmegen93 interaction |36]; dotted line: CD — Bonn interaction [^. Actually the 
Nijmegen93 result is very similar to the Al^lS one and is not reported in this figure, (b) 
The same as in (a), but for B{Q^). (c) The same as in (a), but for T2q{Q'^). Experimental 
data are as in Fig. 1. 



FIG. 4. (a) As in Fig. 3 (a), but for the reduced form factor A{Q'^)/{G']^ ■ F). (b) As 
in (a), but at low Q^. (c) As in Fig. 3 (b), but for the reduced form factor V m^Q"^) I {G"]-)- Fi). 
The Nijmegenl result is very similar to the CD — Bonn one and is not reported in this 
figure. Experimental data are as in Fig. 1 



FIG. 5. (a) The position of the minimum of B{Q^ 
zero of T2o(g^), corresponding to the Gari-Kriimpelmann nucleon form factors 



and (b) the position of the second 

vs the 



nonrelativistic S-state kinetic energy for the deuteron for different realistic interactions. 



FIG. 6. (a) The reduced deuteron form factor A{Q'^)/{G'jj ■ F) obtained with the 
Nijmegen2 interaction for different nucleon form factor models. Solid line: nucleon f.f. of 
Ref. [Q; dashed line: nucleon f.f. of Ref. [^; dotted line: nucleon f.f. of Ref. [^. (b) 
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As in (a), but for the reduced form factor Tm{Q'^)/{G\, ■ Fi). (c) As in (a), but for 72o((5^). 
Experimental data are as in Fig. 1. 
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Fig. la F.M. LEV, E. PACE, G. SALME 
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Fig. lb F.M. LEV, E. PACE, G. SALME 
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Fig. 2a F.M. LEV, E. PACE, G. SALME 
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Fig. 2b F.M. LEV, E. PACE, G. SALME 
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Fig. 3a F.M. LEV, E. PACE, G. SALME 
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Fig. 3b F.M. LEV, E. PACE, G. SALME 
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Fig. 4a F.M. LEV, E. PACE, G. SALME 
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Fig. 4b F.M. LEV, E. PACE, G. SALME 
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Fig. 4c F.M. LEV, E. PACE, G. SALME 
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Fig. 5a F.M. LEV, E. PACE, G. SALME 
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Fig. 5b F.M. LEV, E. PACE, G. SALME 
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Fig. 6a F.M. LEV, E. PACE, G. SALME 
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Fig. 6b F.M. LEV, E. PACE, G. SALME 
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